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ABSTRACT 

A  theoretical   prediction  to  determine  the  frequency  response  of 
a: cascaded  remote  pressure  transducer  system  has  been  developed. 
This  prediction  was  used  to  find  a  theoretical  configuration  which  would 
produce  a  fairly  flat  frequency  response  (1.00  ±  0.05)  over  a  100  hertz 
range.    Such  a  configuration  would  allow  cross -correlation  of  an  unsteady, 
random,  low-speed  flow  field  by  simple  analog  methods. 

The  theoretical  prediction  did  not  hold  true  for  configurations  of 
three  serially  connected  tubes,  but  the  parameters,  namely,  transducer 
volume,  tube  length,  and  tube  diameter,  proved  the  same  trends  experi- 
mentally that  were  predicted  theoretically.     By  observing  these  trends, 
a  solution  of  1.00  ±  0.09  over  a  95   hertz  range  was  obtained. 
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I.      INTRODUCTION 

Because  of  the  increasing  interest  in  measuring  oscillatory  and 
random  unsteady  flow  fields  on  lightweight  aeroelastic  models,  a  need 
has  gradually  developed  over  the  past  several  years  for  a  direct  measuring 
remote  pressure  transducer  system  to  allow  rapid  and  economical  deter- 
mination of  auto-  and  cross -correlation  type  information  plus  the  attend- 
ant spectral  distribution.     Such  a  system  has  been  developed  by  Bergh 
£lj  at  the  National  Aeronautical  and  Astronautical  Research  Institute 
(N.L.R.)  in  Amsterdam,  Holland,  utilizing  precise  lengths  of  tubing  to 
connect  the  remote  pressure  transducers  with  static  orifices  located  on 
the  model  surfaces.     Johnson  \2j  has  confirmed  Bergh' s  one  tube  theory 
and  Allen  l3J  utilized  a  two  tube  configuration:    both  of  these  systems 
required  compensation  for  the  system's  transfer  function  to  be  included 
in  the  data  reduction  process. 

Should  a  configuration  be  found  which  would  produce  a  flat  fre- 
quency response  over  a  specified  range,  cross -correlation  of  an  unsteady, 
random  signal  could  be  achieved  directly  through  the  use  of  an  analog 
circuit.    Since  a  direct  analog  process  does  not  readily  allow  for  cor- 
rection of  the  system's  transfer  function,  the  primary  purpose  of  the 
subject  program  was  to  attain  a  useable  flat  frequency  response  of 
1.0  —  0.05  over  a  0  -  100   hertz  frequency  range.     This  would  allow 
direct  evaluation  of  the  unsteady  aerodynamics  in  many  low-speed  flow 
fields  of  interest. 
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II.      EXPERIMENTAL  PROCEDURE 

A.  GENERAL 

All  experimentation  was  performed  either  in  the  Aero-Mechanics 
Laboratory,  Room  022,  Halligan  Hall,  or  with  the  IBM  360-67  digital 
computer  available  at  the  U.  S.  Naval  Postgraduate  School. 

B.  THEORETICAL  BACKGROUND 

Bergh    [lj    has  developed  a  theoretical  prediction  for  the  pressure 
ratio  of  not  only  a  single  tube  remote  pressure  transducer  system 
(Figure  la  )  but  also  for  a  cascaded  tube  system  employing    N  tubes 
such  as  shown  in  Figure  lb.    The  solution  is  based  on  several  assump- 
tions; the  perturbations  are  of  a  sinusoidal  nature  imposed  on  the  average 
values  of  the  parameters,  fluid  flow  in  the  tubes  is  assumed  laminar, 
and  the  tube  lengths  are,  much  larger  than  their  diameters.     By  introducing 
these  assumptions,  the  governing  equations,  i.e.  ,  the  Navier-Stokes , 
continuity,  energy,  and  state  equations  are  reduced  to  a  set  of  partial 
differential  equations  which  may  be  solved  subject  to  satisfying  the 
boundary  conditions.    The  mass  flow  rates  of  the  air  in  the  system 
between  an/  two  adjacent  tubes  must  be  equivalent.    This  relation  pro- 
vides an  equation  which  reduces  to  the  following  form: 


TV'" 
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where: 


■W^MS?) 


When  a  single  tube  configuration  (see  Figure  la)  is  used,  the  entire 
third  term  of  equation  (1)  disappears,  simplifying  the  equation  to  the 
form  below: 

—   -lco,„($iLi)  +  --^--J-    MfaQj  (2, 

The  subscript  j   refers  to  the  tubes  being  analyzed  in  the  N  tube  system. 
From  equation  (1)    it  is  seen  that  the  pressure  ratio  prediction  is  a  func- 
tion of  tubes  j-1  and  j  +  1   as  well  as  tube  j,  which  in  turn   presents  diffi- 
culties when  optimizing  the  problem  for  the  desired  result.    The  solution 
complexity  is  not  obvious  at  first;  however,  a  closer  inspection  bears 
o.ut  that  all  the  Bessel  functions  are  functions  of  the  imaginary  parameter, 
i  which  are  sometimes  called  Kelvin  functions    \4j    and  hence,  the 

solution  is  a  complex  quantity.    The  complex  pressure  ratio  specifies 
not  only  the  pressure  ratio  magnitude  but  also  the  phase  shift.    The 
equation  also  contains  two  uncertain  quantities;  the  increase  in  trans- 
ducer volume  due  to  diaphragm  deflection,  and  the  polytropic  gas 
constant,    k.     Bergh  proved  both  of  these  terms  to  have  only  minor  effects 
on  the  solutions  to  equations  (1)  and  (2)  and  hence  adiabatic  flow  was 
assumed  while  transducer  volume  changes  were  ignored.    The  system 
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can  be  represented  by  the  acoustic  analogy  provided  by  Olson   [jQ  for 
a  multi-degree  of  freedom  system.    The  analogy  develops  equations  of 
motion  for  a  sound  wave  of  pressure   p  travelling  in  a  tube  and  impinging 
on  a  cavity  volume,  with  the  resistance  represented  by  the  viscosity  and 
the  capacitance  term  dependent  on  the  cavity  volume.     Similarly,  Bergh's 
solution  predicts  the  complex  pressure  ratio  with  the  resistance  and 
capacitance  again  represented  by  fluid  viscosity  and  cavity  volume, 
respectively. 

G.         COMPUTER  TECHNIQUE 

With  the  aid  of  Bergh's  equations,  the  pressure  ratio  at  a  certain 
frequency  can  now  be  predicted.    However,  because  the  pressure  ratio 
of  the  next  higher  tube  is  part  of  the  predicted  ratio  and  the  tube  size 
and  cavity  volume  interact  in  an  irregular  manner,  the  total  solution 
does  not  lend  itself  for  routine  optimizing  techniques  to  produce  a  flat 
frequency  response.    Therefore   the  equations  were  programmed  for  a 
digital  computer  to  facilitate  pressure  response  prediction  using  pro- 
gram-operator chosen  configurations.    The  program  was  designed  for 
use  on  the  IBM  timesharing  system  but  can  also  be  used  for  batch  pro- 
cessing, and  contained  a  program  option  allowing  the  selection  of 
either  a  one-,  two-,  or  three-tube  configuration.     Each  tube  stage  in  a 
multiple  tube  configuration  is  specified  by  three  parameters:  transducer 
volume,  tube  length,  and  tube  diameter.    The  program  was  verified  for 
proper  operation  by  comparing  it  with  both  a  hand-worked  calculation 
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for  a  one-tube  configuration  and  Bergh's    |_1J    results  for  both  one-  and 
two-tube  configurations. 

D.         EQUIPMENT 

1 .         Pressure  Transducers 

The  two  pressure  transducers  were  designed  by 
Professor  L.  V.  Schmidt  [(5J  ,  and  consisted  of  a  Bentley  Detector 
System  mounted  behind  a  0.003-inch  thick  annealed  brass  diaphragm 
which  was  inserted  in  an  aluminum  housing.    The  diaphragm  thickness 
was  selected  to  give  optimum  sensitivity  in  a  i  0.75  psi  range.    The 
detector  system  consisted  of  a  Bentley  316  AJW  proximity  detector 
inserted  into  a  brass  diaphragm  holder.    The  proximity  detector  was  a 
reluctance  gage  which  sensed  change  in  the  air  gap  between  the  reluct- 
ance coil  and  a  metal  target,  in  this  case  the  brass  diaphragm.    The 
diaphragm  deflection  was  directly  proportional  to  an  applied  pressure 
perturbation,  and  the  proximity  detector  output  was  linear  with  respect 
to  the  diaphragm  deflection;  therefore,  the  output  was  a  relative  measure- 
ment of  the  pressure  applied  to  the  diaphragm. 

2  .         Calibration  Chamber 

The  calibration  chamber  was  the  same  as  used  by  Johnson 
j_2]]   and  consisted  of  an  aluminum  block  with  a  1.5-inch  diameter 
cylindrical  chamber  6  inches  long.     The  chamber  ends  were  designed  to 
accept  either  the  transducer  holders  or  the  adapter  plates  which  connected 
the  chamber  to  a  remote  transducer  through  the  transmitting  tubes.     The 
chamber  had  a  pressure  port  to  allow  the  mean  static  pressure  to  be 
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varied  from  an  external  source,  and  also  a  threaded  orifice  designed  to 
accept  either  a  solid  plug  or  a  University  Sound  ID-75  acoustical  driver 
unit.     Both  orifices  were  located  at  the  chamber's  midpoint  so  that 
unsteady  pulses  were  introduced  symmetrically  to  both  transducer 
systems,  namely  the  reference  transducer  at  one  end  and  the  complete 
transducer  system  at  the  other  end. 

E.  CALIBRATION 

1 .         General 

All  calibration  procedures  were  identical  to  those  performed 
by  Johnson    ^2]    and  Allen  \3\  . 

2  .         Proximeter  Calibration 

The  Bentley  316  AJW  proximity  detectors  were  usually  cali- 
brated in  order  to  determine  the  voltage  output  with  respect  to  the  air 
gap  between  the  reluctance  coil  face  and  a  metal  target.     By  mounting 
a  metal  target  on  a  micrometer  stem,  the  target  could  be  backed  off 
known  distances  from  the  reluctance  coil  and  the  output  monitored  on 
a  D.C.  voltmeter.    The  Bentley  D252  carrier  amplifiers  have  a  variable 
amplification  capability  which  may  be  adjusted  until  the  proximity 
detector  output  varies  proportionately  with  the  target  distance.     Because 
of  the  large  apparent  sensitivity  of  the  proximity  detectors  ,  this  pro- 
cedure was  foregone    and  the  detectors  were  calibrated  with  respect  to 
the  pressure  differential  across  the  diaphragm  during  the  static  calibra- 
tion procedures.    The  proximity  detectors  were  inserted  into  the  brass 
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diaphragm  holders  until  their  output  was    -3.0  volts  D.C.:  this  was  the 
range  where  transducer  output  and  diaphragm  deflection  varied  linearly. 

3 .  Static  Calibration 

The  transducer  output  was  connected  to  a  Bentley  D252 
carrier  amplifier  powered  by  a  -18.0  volt  D.C.  regulated  power  supply. 
By  inputing  the  carrier  amplifier  output  simultaneously  with  a  biasing 
voltage  into  an  Astro-Data  D885  differential  amplifier,  the  differential 
amplifier  output,  monitored  on  a  D.C.  digital  voltmeter,  was  adjusted 
to  read  zero  for  zero  pressure  differential  across  the  diaphragm.    A 
U-tube  manometer  and  hand  crank  pump  conveniently  allowed  known 
pressure  differentials  to  be  introduced  into  the  calibration  chamber  for 
static  sensitivity  determination.    The  complete  instrumentation  diagram 
is  shown  in  Figure  2.     Each  proximity  detector  was  calibrated  during 
this  stage  by  introducing  positive  and  negative  pressure  differentials 
of  equal  magnitudes  into  the  calibration  chamber  and  adjusting  each 
carrier  amplifier  until  the  output  was  approximately  the  same  for  equal 
but  opposite  pressure  differentials  over  a    -0.8  to  +0.8  psi  range.     The 
static  calibration  curve  is  shown  by  Figure  3  and  bears  out  the  linear 
sensitivity  variation  as  measured  at  the  carrier  amplifier  output  of  1.37 
and  1.04  volts/psi  for  transducers  1  and  2,  respectively,  over  a  ±0 . 8 
psig  pressure  range. 

4 .  Dynamic  Calibration 

For  the  dynamic  calibration,  a  sinusoidal  perturbation  was 
introduced  into  the  chamber  by  closing  the  static  pressure  port  and 
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inserting  an  ID  75  University  Sound  acoustical  driver  unit  driven  by  a 
Wavetek  114  signal  generator  which  was  amplified  by  a  power  amplifier. 
The  acoustic  driver  unit  was  adjusted  to  operate  at  an  input  level  of 
13.5V  RMS  which  corresponds  to  an  oscillating  pressure  amplitude  in 
the  calibration  chamber  cavity  of  approximately  0.  13  psi  at  120   hertz. 
The  carrier  amplifier  output  was  monitored  with  both  a  true  RMS  meter 
and  a  precision  phasemeter,  and  visually  monitored  on  a  dual  beam 
oscilloscope.     The  dynamic  calibration  instrumentation  diagram  is  shown 
in  Figure  4.     The  frequency  was  varied  from  10  -  12  0   hertz  in  5 -hertz 
increments.     Since  the  disturbance  was  introduced  symmetrically  into 
the  chamber,  there  was  no  appreciable  phase  shift  in  oscillating 
measures  at  either  end  of  the  chamber;  however,  the  sensitivity  of 
transducer  #2  was  0.787  of  transducer  #1  at  every  frequency  increment. 

F.  EXPERIMENTATION 

The  data  runs  were  made  with  the  same  system  as  was  used  for 
the  dynamic  calibration  procedure  with  one  exception;  transducer  #2 
was  remotely  located  from  the  calibration  chamber  through  the  use  of 
nylon  and  stainless  steel  tubing.    The  brass  diaphragm  had  a  centered 
circular  wrinkle  which  allowed  brass  shims  to  be  inserted  between  the 
diaphragm  and  the  plexiglass  spacer  in  the  transducer  housing  without 
hindering  movement  of  the  diaphragm.    The  wrinkle  volume  was  estimated 
as  0.00211  in.      and  the  shims  had  an  inside  diameter  of  0.523  in.: 
this  allowed  any  desired  transducer  volume  to  be  constructed.     Frequen- 
cies of  10  -  120  hertz  were  investigated  in  5-hertz  increments  for 
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different  one-,  two-,  and  three-tube  configurations.    The  ratio  of  trans- 
ducer output  #2  divided  by  transducer  output  #1  was  calculated  and, 
because  the  sensitivity  of  transducer  #2  was  determined  during  dynamic 
calibration  of  transducer  #1  as  0.787,  the  quotient  was  divided  by 
0.787  to  normalize  the  output  of  transducer  #2  with  respect  to  transducer 
#1.     This  determined  the  frequency  response  (transfer  function)  of  a  given 
configuration. 
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III.      DISCUSSION  AND  RESULTS 

The  computerized  solution  of  Bergh's    £l}   prediction  allowed 
quick  and  easy  determination  of  the  frequency  response  for  a  one-,  two- 
or  three-tube  system,  and  could  be  easily  modified  for  a  greater  number 
of  tubes . 

The  one-tube  solution  generally  followed  the  typical  organ  pipe 
solutions  which  produced  sharp  resonance  peaks  at  many  frequencies. 
When  the  resonance  peaks  did  not  occur,  a  heavily  damped  solution 
appeared.    A  solution  was  sought  that  would  reduce  the  resonance  peaks 
or  raise  the  damped  solution  to  the  desired  goal  (1.00  ±  0.05);  however, 
no  such  solution  could  be  found  for  a  0  -  100   hertz  range.    The  data  runs 
did  confirm  Bergh's  theoretical  predictions  for  both  pressure  ratio  and 
phase  shift  for  a  one-tube  configuration,  as  can  be  seen  by  Figure  5. 
Therefore,  a  two-tube  solution  was  examined  in  hopes  that  the  pressure 
of  an  additional  tube  would  counter  the  undesirable  characteristics  of 
the  first  tube.    Such  a  configuration  was  found  with  the  desired  frequency 
response  of  1.00  ±  0.05  from  10  to  100  hertz  as  shown  by  Figure  6. 
However,  it  was  felt  that   the  first  tube  in  the  system  had  a  diameter 
that  would  result  in  the  static  orifice  on  the  body  being  large  enough  to 
perturb  the  flow  in  the  region  of  static  pressure  measurement.    The  two- 
tube  configuration  also  verified  Bergh's  theoretical  prediction  of  phase 
shift,  but  the  pressure  ratio  was  consistently  lower  (6%)  than  predicted 
for  frequencies  greater  than  40  hertz  (see  Figure  7).    A  three-tube 
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configuration  was  then  sought  in  order  to  reduce  the  diameter  of  the  first 
tube  to  an  acceptable  value  (less  than  0.05  inch).    A  solution  which  met 
all  desired  specifications  (1.00  ±  0.03  over  a  10  -  100  hertz  range)  was 
found  (see  Figure  61!  and  physically  constructed  in  order  to  obtain  experi- 
mental results  to  confirm  the  configuration  selected.    The  configuration 

specifications  were:  _   .  m   ,  m  , 

Tube  Tube  Tube 

#  1  #  2  #3 

Transducer  volume  (in.  3)  o.O  0.0  0.0238 

Tube  length  (in.)  0.45  29.0  21.0 

Tube  diameter  (in.)  0.0465  0.078  0.03 

The  three-tube  solutions  again  confirmed  the  phase  shift  prediction; 
the  experimental  values  were  within  7%  of  their  predicted  values  (see 
Figures  8a-d).       However,  as  was  the  case  with  the  two-tube  solution, 
the  three-tube  experimental  configuration  produced  pressure  ratios  less 
than  those  predicted  by  Bergh's  solution  for  frequencies  above  30  hertz 
with  the  values  lower  by  as  much  as  30%  (see  Figures  8a-d) .     This 
characteristic  made  the  theoretically  acceptable  solution  (Figure  8a)  of 
limited  practical  value  since  the  experimental  response  was  consider- 
ably below  the  prediction.     Even  though  the  actual  solutions  were  not 
as  predicted,  Figure  10  shows  trends  comparable  to  those  predicted  by 
Bergh's  solutions  for  identical  configurations  (Figure  9). 

From  Figures  8b-d,  it  was  observed  that  changing  the  lengths  of 
tubing  raised  or  lowered  the  overall  curve  and  a  decrease  in  transducer 
volume  (compare  Figures  8a,  b)  tended  to  increase  or  raise  the  overall 
curve  in  the  region  of  interest  and  decrease  the  slope  above  80  hertz. 
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A  small  volume  was  tested  and  a  possible  solution  of  1.00  ±  0.09  over 
a  10-95  hertz  range  was  finally  attained  (see  Figure  8d) .     This  configu- 
ration appears  acceptable  for  direct  analog  data  processing  in  low-speed 
fields  of  interest.     Much  of  the  separated  flow  phenomenon  is  involved 
at  non-dimensional  frequencies  or  Strouhal  numbers,  S,  less  than  0.4 
which  would  typically  correspond  to  frequencies  less  than  100  hertz  for 
a  five-inch  diameter  circular  cylinder  at  a  tunnel  speed  of  100  feet  per 
second.     The  solution  configuration  is  specified  below: 


Tube 

Tube 

Tube 

#1 

#  2 

#  3 

0.0 

0.0 

0.0174 

0.45 

28.5 

21.0 

0.0465 

0.078 

0.03 

Transducer  volume  (in.  3) 
Tube  length  (in.) 
Tube  diameter  (in.) 

In  addition  to  the  tube  lengths  and  transducer  volumes,  the  three 
tube  diameters  were  very  important  parameters  when  attaining  the  theo- 
retical response  desired.    The  diameter  effects  and  trends  were  not 
examined,  however,  due  to  time  and  facility  limitations  and  particularly 
due  to  the  limited  selection  of  commercial  tube  diameters  available  with 
acceptable  inside  diameter  specifications. 

The  purpose  for  attaining  a  flat  frequency  response  was  to  allow 
easy  auto-  and  cross-correlation  of  a  body  in  a  low-speed  unsteady 
flow  employing  a  simple,   low  cost  analog  method  shown  by  Figure  11. 
The  cross-correlation  coefficients,  which  relate  the  coherence  of  two 
physical  measurements,  say  at  stations  (1)  and  (2),  are  in  general 
defined  for  a  specific  value  of  time  lag,     L  ,  by: 
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However,  the  elementary  analog  correlation  circuit  shown  in 
Figure  11  does  not  have  a  time  lag  capability  and  hence  the  benefit  of 
simplicity  is  associated  with  a  loss  of  information.    Therefore,  the 
information  to  be  obtained  must  be  at  zero  time  lag  as  specified  by 
equation  (3): 


Since  the  system  in  consideration,  i.e.  ,  a  remote  pressure  trans- 
ducer system,  can  be  likened  to  a  typical  transfer  function  system  shown 
by  Figure  la,  b,  the  cross-correlation  inputs,  e.    and  e    ,  can  be  repre- 
sented by  the  transducer  outputs ,  p     and  p    .     If  the  disturbance  were 
of  a  sinusoidal  nature  as  it  was  for  the  laboratory  experiment,  then  the 
two  measuring  stations  become  strictly  functions  of  the  frequency, 

Suppose: 

-/°i  =  -po\     sin     co  t 
~P^~  -p*-L   Sin    (cot-©) 

-£.  (+)  =  i°oi  Gr^cS)    sir)   (ojf-c/),) 


then: 


37 


Substitution  into  equation    (3)   will  yield  a  cross -correlation  coefficient 
as  follows: 

d    (o)  -  G>  (c°^  G^  ^  cos  ^ e  4  ^  z  <h  ) 

From  equation  (4)    it  can  be  seen  that  for  the  case  of  a  pair  of  information 
signals  at  a  unique  circular  frequency,  co,  with  the  respective  transfer 
functions  G     (co)  and  G     (co)  equivalent,  then: 


Rit.  (°)  -   COS    © 


and  a  flat  frequency  response  is  not  necessary.     It  is  only  necessary  to 
have  identical  transfer  functions  used  with  both  signals. 

However,  this  analysis  eliminates  one  of  the  most  important  types 
of  flow  fields  ,  the  unsteady  random   signal  such  as  that  produced  by  flow 
separation  from  a  bluff  body.     In  such  cases  the  cross -correlation 
coefficient  can  no  longer  be  simplified  to  equation  (5).    Instead,  because 
of  its  reliance  on  G  as  a  function  of  tO  ,    the  output  of  the  cross- 
correlation  system  becomes: 

(6) 


fU*)  = 


where 

<P         =  power  spectral  density  of  signal  1 

T  22   =   P°wer  spectral  density  of  signal  2 

co-spectral  density  of  signals    1  and   2 


f 
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At  the  input  of  the  transducer  system,  the  cross -correlation  coefficient 
is  not  a  function  of  G  and  thus  reduces  to  the  form: 


R    (o)_         5*    1^0^ 


In  order  for  (6)  to  equal  (7)  for  any  random  signal,  the  following  relation- 
ship is  obviously  necessary. 

Gr  fo)    =  G2  (w)    =    l.O  (8) 

Should  (8)  be  the  case,  then, 

R..(0)   .  =     R      (0)        .      . 

12        input  12         output 

and  the  system  ho  longer  relies  on  the  gain  of  the  transducer  system. 

A  complete  derivation  of  the  cross-correlation  equations  may  be  found 

in  Appendix  C . 

Because  of  the  flat  frequency  response,  the  collected  data  may  be 

inputted  directly  to  the  cross -correlation  circuit  shown  in  Figure  11  and 

the  output  monitored  by  either  a  D.C.  voltmeter  or  an  x-y  plotter  to 

I  i 

0        2  9        2 

obtain  quantities      <e  e   >    ,  ^e.  )    ,    and  \e    }    ,  which  allows  a  simple 
calculation  of  the  cross-correlation  coefficient  at  zero  time  lag. 
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IV.     CONCLUSION  AND  RECOMMENDATIONS 

Although  the  breakdown  of  Bergh's  theory  seems  to  be  directly 
proportional  to  the  number  of  tubes  in  the  configuration,  it  does  predict 
fairly  accurately  the  trends  that  each  parameter  impresses  on  a  system. 
There  is,  therefore,  a  good  possibility  of  obtaining  a  flatter  frequency 
response  than  the  one  obtained  in  Figure  8d.    Time  limitations  required 
the  termination  of  the  experiment  and  the  three-tube  configuration  was  in 
no  way  completely  examined. 

It  should  be  noted  that  the  flat  frequency  response  traits  considered 
herein  applied  to  a  fairly  restrictive  range  of  frequencies,  namely  less 
than  100  hertz.    The  processing  of  a  broader  class  of  signals,  i.e.  ,  a 
higher  range  of  frequencies,  will  probably  require  the  utilization  of  several 
steps  or  procedures,  including: 

a.  Recording  of  the  analog  form  of  the  data  signals  on  a  magnetic 
tape  recorder. 

b.  Digitizing  the  analog  signals  using  a  common  time  base  or 
clock,  and 

c.  Processing  the  digitized  data  by  means  of  a  digital  computer. 
Portions  of  this  procedure  have  been  within  the  scope  of  the  school's 

experimental  facilities  within  recent  times,  cf.  ,  Jesberg  [yQ  ;  however, 
new  equipment  additions  or  replacements  make  this  former  technique  no 
longer  practical  until  a  repeat  is  done  of  the  work  by  Jesberg  \7J  . 
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It  is  observed  from  Figures  7  and  8a-d  that  the  theory  under- 
predicts  the  pressure  ratio  for  both  two-  and  three-tube  configurations, 
yet  the  phase  prediction  is  very  close  to  the  correct  value. 
The  following  recommendations  are  therefore  made: 
1.)       Further  experimentation  of  parameter  effects  in  order  to 

obtain  a  flatter  frequency  response. 
2.)       Possible  modification  of  the  computer  program  to  allow 
compensation  for  the  difference  between  the  actual  and 
theoretical  response. 
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APPENDIX  A 

THEORETICAL  RESPONSE  OF  A 
REMOTELY  LOCATED  PRESSURE  TRANSDUCER  SYSTEM 

Bergh's  theoretical  model  applies  to  any  number  of  transmitting 
tubes  in  a  cascaded  tube  system  as  shown  by  Figure  lb;  however,  many 
assumptions  are  made  in  order  to  reduce  the  governing  equations  to  a 
workable  form.     These  assumptions  are  as  follows: 

a)  the  sinusoidal  disturbances  are  very  small: 

b)  the  inner  diameter  of  the  tubing  is  small  with 
respect  to  its  length: 

c)  laminar  flow  exists  throughout  the  system: 

d)  the  unsteady  disturbances  are  perturbations 
imposed  on  the  average  values  on  the  parameters 
indicated  below: 

Pressure 

Density 

Temperature  (Al) 

Axial  Velocity 

Radial  Velocity 

These  assumptions  allow  the  governing  equations  to  be  reduced  to  the 
following  forms: 


1) 

P 

= 

Po 

+ 

iwt 
pe 

2) 

P 

= 

Po 

t 

iwt 
pe 

3) 

f 

= 

T 
o 

JL 

m  iwt 
Te 

4) 

u 

= 

u 
o 

+ 

iwt 
ue 

5) 

V 

= 

V 

o 

+■ 

iwt 
ve 
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Navier-Stokes  Equations 
Continuity  Equation 


O  (A2a) 

2>v 


(A3) 


Perfect  Gas  Law 


■5*  t '  +  ^  -1  =  f 


(A4) 


Energy  Equation 


(A5) 


Since  the  equations  are  now  partial  differential  equations,  boundary 
conditions  can  now  be  specified: 

a)         At    r  =  R,  the  tube's  inside  wall, 

1)  u(X,R)    =   0 

2)  v(X,R)    =   0 

3)  T(X,R)    =   0 
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b)         At  r  =  0,   the  tube  center, 

1)  v(X,0)   =  0 

2)  u(X,0)   =  finite 

3)  T(X,0)    =   finite 

4)  p(X,0)    =   finite 

5)  p(X,0)    =   finite 

From  equation  (A2a)    it  is  evident  that   p  =   p(x),    which  allows  a 
separation  of  variables  technique. 

assume: 

T=  f(x)  K(2) 

where: 


(A6) 


Vi. «  /  p*^y/z. 


/x_  y 


A 


R 


2=  T  p' 

Equation  (A5)  takes  on  a  Bessel  equation  form  of  the  following  type  in 
this  type  of  solution: 

■p       i 


M*)a  c,  •£  C«)  +  <^  X  (*)  +■ 


(OfcC,.   f(x)  (A7) 


Apply  boundary  conditions  a3  and  b3 

I  -PC*) 


f  (*)  =  - 


c,o;^p^)    po^cp 


(A8) 
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Using  (A7)    and    (A8)  in  (A6)  yields: 


"SF* 


I" 


>r-i 


I- 


(A9) 


Similar  methods  using  boundary  conditions  al   and  b2  can  be  applied  to 
solve  equation  (1)  as  given  below: 


Us 


3-.(£«) 


-I ,- 


(A10) 


Substituting  equations  (A9^  and  (A10),   integrating  with  respect  to  r, 
and  applying  boundary  condition  a2  will  put  equation  (A7)  in  the  form: 

_£   <d*P     JV.(«0 
2-     dx^    T0K) 


Since  v(X,0)    =   0   implies  F(X)    =   0, 


CO 


'-X-p 


1+ 


*  j;  («•  ft-  "*) , 


=  o 


The  above  equation  can  be  solved  for  the  homogeneous  solution: 

-p=  Ae**  +  Be"** 

where  ,  y  1/ 


°-»  Ji'M«) 


*  Jo    («PA) 
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The  coefficients  A  and   B  can  be  determined  only  after  the  boundary 
conditions  for  both  ends  of  the  tube  have  been  specified.     It  should 

also  be  noted  here  that  the  Bessel  functions  in  the  equations  are  complex 

3/o 

numbers  since      o^  =  o(,  ( i  y).    This  is  discussed  further  in  Appendix  B 

Two  additional  assumptions  must  now  be  made  if  the  solution  is  to 
cover  an  N-tube,  N  volume  configuration: 

c)  '  the  pressure  and  density  in  the  instrument  volumes 

are  time  dependent  only 

d)  the  pressure  expansion  in  the  instrument  volume 
is  a  polytropic  process  described  by 


=    constant 


*K 


The  constants  A  and  B  are  now  determined  in  the  following  procedure 
For  any  tube  j , 


at   x    =  0 

J.  (All) 


p.    ,    =  A.   +  B. 

J-l  J  J 


at    x.  =   1. 
J  J 


-ft  =  (l;C^    +    Bie"**L* 


)    ["n    ~foi      „    --foil 

(A12) 


•■^{^^^[v^^] 
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Equations    (All  and    (A12)  allow  solution  for  A.   and   B.. 


6,  =  &"  € 


3      e+^>-e-^i 


The  mass  leaving  the  tube  can  now  be  calculated: 


/r> 


,.=  £ 


uz  zvrrdv  = 


X  to  3"0(c*v  ) 


fl^e^-B^e-***' 


For  any  tube   j  +    1 , 


at   x.   _   =  0 


p.  =  A.    .    +  B.    _ 

J  J  4-1  )  +  l 


£<iii^a+i 


The  mass  leaving  tube   j  +  1    is  now  calculated: 


>t  do 


^"o  (<*$+») 


(A13) 


(A14) 
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Equations    (A13)    and    (A.14)  allow  solution  for  A.    n    and   B.    , 

j+1  j+1 


*i*r 


v^-nc**'1*' 


e 


4vm£.*h  -  e"<'>>+,L>+' 


The  mass  variation  within  the  instrument  volume  must  equal  the  differ- 
ence in  the  mass  leaving  tube   j   and  entering   j+1,    namely: 


cJ  WV 


:*,+ 


-  =  c^V^*  Vvi  (°5+kOs  ^•"^•) eXaj       (A15) 

■J"  '    O        Q»6  » 


Substitution  of  the  pressure  coefficients  into    (A15)   and  recalling  the 
definitions 

cosh  4>t     =   and   sinh  cot    -    

the  final  pressure  ratio  prediction  is  obtained. 

+  ^^i-«-»     fta+i   Ek^ili)  ^^^     S/NH  (j^jj    x 


{ 


««(*jtlli«)-^ 
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APPENDIX  B 

COMPUTERIZATION  OF  THE  THEORETICAL  RESPONSE 
OF  A  REMOTE  PRESSURE  TRANSDUCER  SYSTEM 


A.         Computer  Program  Development 

In  Bergh's  theoretical  development  (Appendix  A) ,  the  Bessel 

3/2 
function  parameter,     <*,  contains  the  imaginary  term,  i  '     ,  which  can 

easily  be  represented  by  the  Kelvin  functions  Ber,  Bei,  Ber   ,  and  Bei    . 

These  functions  were  used  in  the  program  in  the  form  of  an  infinite 

power  series  and  terminated  when  the  terms  become  smaller  than  10 

The  series  expansions  were  taken  from  McLachlan     \4 J     and  altered  into 

the  form  below: 


J~o6*)  -    ber  y    +  X    bei    y 

//°o6j\'/2. 

[jz)    =  ,*.'*  y 


(Bl) 


«*j*R<*£r*** 


ber    y    =   )-  -~-  -h— ?— 


■f... 


l  Yx         Yc  Y' 

bei    V  =  — - —    -i 


49 


Cf,  (<*)  =    ber,   y    +  jl     bel  ,  y 

y    r         y"2-  yf  y6 

k  /  r     yz        y*  y  * 


(B2) 


y 


(B3) 


B.         Computer  Program  Operation 

The  computer  program  was  designed  primarily  for  use  on  a  time- 
sharing system  remote  terminal  where  quick  output  could  be  obtained 
for  a  given  input.    However,  it  may  also  be  modified  for  batch  processing 
as  explained  below.    The  computer  program  allows  selection  of  either 
a  1,  2   or  3  tube  configuration  by  the  operator  and  prints  out  the  pressure 
ratio  and  phase  shift  for  multiples  of  5   hertz  from  10  to  120  hertz.    The 
flowchart  is  given  in  Figure  (Bl). 
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/  READ     \ 

/v3,   t3,  dx 

V2.     T2,    D2 
Vl'     Tl<    Dl 

\ 

' 

#3    TUBE 
.    ANALYSIS 

1,2,   OR   3   TUBE 
CONFIGURATION    ? 


V. 


READ 

V    D2 
Ti  ,    D, 


# 2    TUBE 
ANALYSIS 


#  1    TUBE 
ANALYSIS 


WRITE 
PRESSURE   RATIO 
&  PHASE   SHIFT 


READ 

v  .  T  .  Dj 


yes 


ADDITIONAL^ 
SOLUTION 
DESIRED   ?. 


FIGURE   Bl 


COMPUTER  PROGRAM   FLOWCHART 
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APPENDIX  C 
CROSS -CORRELATION  OF  A  RANDOM  UNSTEADY  SIGNAL 

The  cross -correlation  coefficient  may  be  defined  in  its  most 
general  form  similar  to  methods  described  in  Bendat  and  Piersol   £8"}. 

R   \+)  =  <e,(+) e^+-^r)> 

However,  because  of  the  zero  time  lag  limitation  of  the  analog  circuit 
shown  in  Figure  11,  this  reduces  to  the  form: 


R„.fe),         <e'^ 


(C2) 


For  the  system  input,  the  following  terms  are  defined  as: 

•*0       ili        '  (C3) 

<e<2  >  =    57°  $  „  (to)  clco 
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d  =  power  spectral  density  of  signal  1 
d>  _  =  power  spectral  density  of  signal  2 
<£-.  9   =  co   spectral  density  of  signals    1   and   2 

After  passing  through  the  system,  the  quantities  have  a  transfer  function 
imposed  on  their  true  output,  and  the  following  terms  are  defined: 

e.,(">)  =  A(e-^-'aj+  <C4) 


thus: 

w>  =  ir^M  j,Hd«  (C5) 


A  relationship  between  the  output  and  input  cross -correlations  may  be 

obtained  by  direct  substitution  of  (C3)  and  (C5)  into    (C2),  yielding: 
©o      ^ 
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There  are  two  possible  simplifications  available  at  this  point.     If  the 
input  were  of  a  harmonic  nature,  say  ac  frequency     Co  ,  then  the  expres- 
sions for  spectral  density  would  simplify  and  yield: 

(C7) 

^e^  =  C°  gjlfr»)S(«'-<Ui)'<'&j  =  5ltCoj,) 

similarly: 

(C8) 

and  the  dependence  of  the  zero  time  lag  cross -correlation  coefficient 
upon  the  system's  transfer  function  would  disappear  due  to  cancellation. 

However,  if  the  input  were  of  a  random  nature,  then  the  only 
simplification  readily  available  is  to  have: 
G     (to)    =  G2  (cu)    =    1.0 

in  which  case  the  original  input  cross -correlation  is  obtained  since  the 
system  transfer  function  factors  out  from  the  integrands,  cancels,  and 
simplifies  to  equations  .* 

<$,'•>=   <e,l>  =51°  $„(«-)    dw 

<£*  >  =  <e*>  =  ST  £„.«<»)  dw  (eg) 
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<a,n>«<«.«o  »5Tflv(^a 


Thus: 

£/z.   Co)    -     #^  (o)  (CIO) 

Since  the  latter  case  also  applies  to  the  harmonic  signal  data 
processing,  this  is  a  most  desirable  situation  and  is  the    motivation 
behind  seeking  a  flat  frequency  of    1.00  ±  0.05. 
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